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ABSTRACT 

It has become increasingly apparent that traditional hydrodynamical simulations of 
galaxy clusters are unable to reproduce the observed properties of galaxy clusters, 
in particular overpredicting the mass corresponding to a given cluster temperature. 
Such overestimation may lead to systematic errors in results using galaxy clusters as 
cosmological probes, such as constraints on the density perturbation normalization erg. 
In this paper we demonstrate that inclusion of additional gas physics, namely radiative 
cooling and a possible preheating of gas prior to cluster formation, is able to bring 
the temperature-mass relation in the innermost parts of clusters into good agreement 
with recent determinations by Allen, Schmidt & Fabian using Chandra data. 
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1 INTRODUCTION 


clusters. For example, the recent Chandra results of Allen 


Reprod ucing the observed number density of rich galaxy 
clusters has long been thought to be one of the most re¬ 
liable constraints on the matter power spectrum on short 


Schmidt & Fabian (2001, hereafter ASFOl) indicate that, at 


(Evrard 198£; Henry & Arnaud 1991 

; White, Efstathiou & 

Prenk 1993|; Eke, Cole & Frenk 1996 

fViana & Liddle 19961 

1999; Henry 1997, 2000; Blanchard et al. 2000; pierpaoli. 


ically yielding ns ~ 0.9 to 1.0 for the currently-favoured 
ACDM model with matter density flo — 0.3. However, re¬ 
cently evidence has begun to accumulate from a number of 
sources that this may be a significant overestimate, perhaps 
by tens of percent. For exam ple, the required ag estimated 


from the 2dF galaxy survey (Lahav et al. 2002; Verde et al 


2002), or from that survey combined with other probes (Efs- 


tathiou et al. 2002), is significantly lower, and there are now 


suits (Borgani et al. 2001; 

Reiprich & Bohringer 2002 

; Viana, 

Nichol & Liddle 2002). 


who used an ob served relation between cluster tempera - 


ture and mass (Finoguenov, Reiprich & Bohringer 2001) 


rather than one derived from hydrodynamical simulations. 
This last result is particularly significant, and points to the 
increasingly evident result that traditional hydrodynamical 
simulations, which include only adiabatic gas heating during 
collapse, are unable to reproduce the observed properties of 
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least in the inner regions where data exists, clusters are con¬ 
siderably hotter for a given mass than predicted by adiabatic 
simulations. 

Here we address the question of whether the inclusion of 
additional gas physics, both radiative cooling of the gas and 
preheating of the gas before cluster formation, is capable of 
bringing the simulations into agreement with observations. 
We concentrate only on the inner regions of clusters, for 
which temperature profiles have been measured by Chan¬ 
dra, and we find that indeed the observations can be repro¬ 
duced. In itself this is not sufficient aid to theorists seeking 
to constrain as, which requires an accurate description of 
clusters out to the virial radius, but this encouraging result 
suggests that simulations may soon be useful for this pur¬ 
pose. We will explore the cluster temperature-mass relation 
out to larger radii in a forthcoming paper. 


2 DETERMINING CLUSTER MASSES 

2.1 Using the Hydrostatic Equation to measure 
masses 

The distribution of hot gas in a cluster can be used to mea¬ 
sure its mass. The intracluster medium is generally assumed 
to be in hydrostatic equilibrium in a spherically-symmetric, 
static potential and the mass is determined by balancing 
pressure support against the gravitational attraction: 
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( 2 ) 

is the pressure of the intr- 


aclustei l medium. Pgaa its density, k is the Boltzmann con- 


stant, G the gravitational constant, and /rmn — 1.0 x 10“ g 
the mean mass per particle. 

Where Pgas is simply the thermal pressure, then Pgas 
corresponds to the gas temperature. One should really in¬ 
clude all forms of pressure support for the gas: kinetic 
(i.e. turbulence), magnetic, coupled relativistic particles, 
etc.. We have checked in our simulations that the contribu¬ 
tion from kinetic motions within the inner regions of clusters 
is small and so for clarity of presentation we stick with the 
thermal pressure in this paper. 

Note that the determination of the mass within radius 
r depends only upon the properties of the gas at that ra¬ 
dius; in particular the mass determination is not affected by 
conditions in the outer parts of clusters whe re th e properties 
are poorly determined. We show in Section 4.1 that the use 
of equation (^) leads to good estimates of the mass within 
regions accessible to Chandra observations. This lends cre¬ 
dence to the mass determinations from X-ray observations 
that we use in this paper. 


2.2 Observed and simulated temperature—mass 
relations 

If we apply equation to measure the mass, Ma, within 
a radius, va, for which the enclosed density is A times the 
critical density, then 


kT{rA) 

keV 
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where Va = — din Pgas/din r measured at ta- 

For self-similar clusters, for which Va is independent 
of mass and T^a) is a constant multiple of the observed 
temperature, Tx, we recover a scaling relation Tx oc Ma^^- 
In order to compare different results we will write 


kTx 

keV 


= ^A 


Ma 
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(4) 


Aa gives the relative normalization of the relations; if 
a 7 ^ 2/3 then this comparison will be exact only for 
Ma = lO^®/i“^M 0 . Theorists can most easily predict the 
mass within the virial radius of cluster, corresponding to 
A ~ 111 for our chosen cosmology, but X-ray observations 
do not extend to such large radii and so some degree of ex¬ 
trapolation (i.e. the extension of a mass-model beyond the 
range of the observational data) is usually required. 

A useful su mmary of simulated and observational re¬ 
sults is given in Afshordi fc Cen (200^). The main simu¬ 
lated cataloguesaxe^^^^^^^^!"7l9^)) Thomas et al 


(2001) and Bryan & Norman (199f). The former two, us¬ 


ing smoothed-particle hydrodynamics (SPH), found A 200 ~ 
8.0, whereas the latter, using an Eulerian grid-code found 
A 200 « 7.0. In a more recent paper, in which they consid¬ 
ered a variety of definitions of X-ray temperature for an 


ensemble of 24 highly-resolved SPH clusters, Mathiesen & 


Evrard (2001) agreed with this lower normalization. None 


of these simulations included radiative cooling. 

By contrast the observational results using resolved 
surface-brightness and temperature p rofiles from ROSAT 
and ASCA have higher normalizations. Xrq^jIiii^^Vu_^2001 ) 
find A 200 ~ 9.4, whereas Horner , Mushotzky & Schari 
(1999|) and Finoguenov et al. (2001[) both find A 200 ~ 11.2. 


Despite the uncertainty in both the simulated and the 
observational results, it is clear that the observed normaliza¬ 
tion of the temperature-mass relation significantly exceeds 
the simulated one. 


2.3 Observational determination of M 2500 

In this paper, we attempt a partial resolution of the differ¬ 
ences between simulations and observations discussed in the 
previous section. In doing this, we concentrate on the results 
of ASFOl. The reasons for this are fourfold: 

(i) Chandra observations give the best available estimates 
of density and temperature profiles for the X-ray emitting 
gas. 

(ii) The mass estimates are, mostly, backed up by obser¬ 
vations and modelling of gravitationally-lensed background 
galaxies (arcs). 

(iii) They present results for the mass-weighted temper¬ 
ature of the gas. The use of mass-weighted rather than 
emission-weighted temperatures greatly simplifies the com¬ 
parison of simulations and observations. 

(iv) They do not attempt to extrapolate their results be¬ 
yond the radius that is accessible to observations. 

ASFOl measured the mass and temperature of 5 clusters 
within r 25 oo. They found a best-fit slope for the mass- 
temperature relation that is consistent with the self-similar 
value of 1.5. We rewrite their relation here with mass as the 
ordinate as we are complete in mass rather than tempera¬ 
ture in our simulations: 
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where T 2500 is the mass-weighted gas temperature within 

’' 2500 - 


3 SIMULATIONS 

The simulations that we discuss in this paper were carried 
out using the Hydra A-body/ h ydrodynamics code (Couch- 
man, Thomas & Pearce 1995; Pearce & Couchman 1997) 


on the Cray T3E computer at the Edinburgh Parallel Com¬ 
puting Centre as part of the Virgo Consortium programme 
of investigations into structure formation in the Universe. 
They will be described fully in a longer paper (in prepara¬ 
tion) and so we just summarize the properties here. Note 
that the simulati ons are very similar to those discussed in 
an earlier paper, Muanwong et al. (2001), but the parame¬ 


ters have been slightly adjusted so as to give a better fit to 
the observed luminosity-temperature relation of clusters. 

In this paper we present results for a single cosmol¬ 
ogy with density parameter flo = 0.35, cosmological con¬ 
stant Hao = 0.65, baryon density flbo = 0.038, Hubble 
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Figure 1. The solid line shows the actual mass distribution in 
the cluster; the dashed line shows the mass from equation 
using the thermal pressure. 



'’/'’2500 

Figure 2. Temperature profile of the gas within r25oo for on® 
example cluster in the three simulations. 


parameter h = 0.71, cold dark matter density fluctuation 
parameter T = 0.21, normalization as — 0.90 and gravita¬ 
tional softening s = 25/i“^kpc. 160® particles each of gas 
and dark matter were used in a box of side 100h“^Mpc, 
giving particle masses of mgas ~ 2.6 x lO®/i“^M0 and 
m-dark ~ 2.1 X lO^®h~^M0, respectively. 

Three simulations were undertaken, each with identical 
initial conditions so that the clusters that they produce are 
directly comparable, but with different cooling properties: 


Non-radiative: An adiabatic run includ ed for testing and 


comparison purposes only. As discussed in Muanwong et al 


(2001^the gas in the centres of the clusters extracted from 
this run has short cooling times and would not be present in 
real systems. This simulation vastly overestimates the X-ray 
luminosity of the clusters. 

Radiative: Th is run includes radiative coo ling using the 
cooling tables of Sutherland & Dopita (1993) and assumes 
a time-varying metallicity Z = 0.3 {t/to)ZQ, where t/to is 
the age of the Universe in units of the current time. Cooled 
material is permitted to form stars, removing low-entropy 
material with short cooling times from the centres of the 
clusters. As shown in Pearce et al. (200C), the inflowing gas 
which replaces it has a higher net entropy and hence a higher 
temperature. A drawback of this run is that the cooling is 
limited only by numerical resolution and is therefore very 
ad-hoc. In this run up to half the baryonic mass in clusters 
has cooled, much more than is observed. 

Preheating: This run also includes radiative cooling. Ad¬ 
ditionally, at z = 4 the energy of each gas particle was in¬ 
creased by an amount kT = 1 keV to crudely model energy 
injection at high redshift. This has the effect of strongly sup¬ 
pressing cooling so that by the end of the simulation under 
1 per cent of the gas in clusters has cooled, much less than 
is observed. 


The Radiative and Preheating runs both reproduce the ob¬ 
served luminosity-temperature relation whilst having very 
different amounts of cooled gas. One might hope, therefore, 
that their thermal properties also bracket those of real clus¬ 
ters. We justify this statement further in Section 4.3 below. 


4 RESULTS 

4.1 Hydrostatic equilibrium 

Figure compares measures of enclosed mass versus radius 
for the third-largest cluster in the Preheating simulation 
(similar results are obtained for the Non-radiative and Ra¬ 
diative simulations). This particular cluster was chosen be¬ 
cause the two largest clusters both show signs of disturbance 
within rsoo. The solid line shows enclosed mass as a function 
of radius. The dashed line shows the mass estimated from 
equation (^) using the thermal pressure; a similar result is 
obtained using the total (thermal plus kinetic) pressure but 
we have omitted this from the plot for clarity. 

The estimated mass jiggles up and down because of vari¬ 
ations in the local pressure gradient. This scatter could be 
much reduced by the fitting of a smooth curve to the pres¬ 
sure profile, but for the purposes of this paper we merely 
wish to make the point that the clusters are in approxi¬ 
mate hydrostatic equilibrium within r25oo- This holds true 
also for the largest two clusters and for all others that we 
have tested. Note, however, that many of these other clusters 
show significant departures from equilibrium within r2oo- 

We conclude that the use of the Equation of Hydrostatic 
Support to measure masses from X-ray observations of the 
intracluster medium within r25oo is likely to be accurate to 
within about 10 per cent with no systematic offset to high 
or low masses. 


4.2 Temperature profiles 


Figure g shows the projected temperature profiles within 
r2500 of the third-largest cluster in each of the three sim¬ 
ulations. Note how the inflow of higher entropy material 
has raised the temperature of the gas within r25oo in both 
the Radiative and Preheating simulations relative to that in 
the Non-rad iative simulation, an effect first noted by Pearce 
et al. (200(]|) . The temperature profiles in this Figure differ 
from those shown in Figure 1 of ASFOl in that they decline 
slightly beyond 0.25 r25oo rather than remaining isothermal. 
There are other clusters in our sample, however, for which 
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Figure 3. Entropy profile of the gas for one example cluster in 
each of the three simulations. 


the temperature profiles are very similar to those seen in the 
observations. 


4.3 Entropy profiles 


Figure 4. This panel shows the mass-weighted temperature- 
mass relation for gas within r 2500 for clusters extracted from the 
Non-radiative simulation. Different symbols correspond to clus¬ 
ters with different amounts of substructure, as discussed in the 
text. The solid line shows the best-fit from ASFOl; t he dashed line 
shows fextrap olated! results from the simulations of Mathiesen & 


Evrard (2001). 


The differences between the three simulations are best de¬ 
scribed in terms of their entropy profiles. As Figure^ shows, 
the entropy within riooo in the third-largest cluster is raised 
in the Radiative and Preheating runs relative to that in the 
Non-radiative run. This results in a higher gas temperature 
and a less-peaked density profile. 

Although the Radiative and Preheating runs have very 
different cooled mass fractions, it is not surprising that they 
have similar entropy profiles as we have adjusted the numer¬ 
ical resolution and feedback parameters so as to reproduce 
the observed X-ray luminosity-temperature relation and it 
has been recognized for some time that this requires an 


excess of en t ropy in clus ter cores (Evrard & Henry 1991 
Kaiser 1991; Bower 1997| ). Once the entropy profile is fixed 
then, assuming hydrostatic equilibrium, the cluster temper¬ 
ature profile is uniquely determined. As it happens, both 
the Radiative and Preheating runs give very similar entropy 
profiles and hence very similar temperature profiles for this 
cluster (and similar results are obtained for other clusters). 
This gives us confidence that other models that reproduce 
the luminosity-temperature relation, and in particular ones 
that give the correct cooled gas fraction, would have similar 
thermal properties to those discussed here. 


4.4 Temperature—mass relation 


The temperature-mass relation for the clusters extracted 
from the Non-radiative simulation is shown in Figure [f|. The 


dashed line is the relation from the simulations of Mathiesen 


&Z. Evrard (2001), extrapolated from their results for an over¬ 


density of 500 as described in ASFOl. Given this extrapola¬ 
tion our results are in good agreement with theirs over the 
range for which our temperatures coincide, kT > 1.5 keV. 
By contrast, the relation for observed clusters obtained in 
ASFOl, shown by the solid line, has a significantly lower 
mass normalization for a given temperature. 


Figure ^ shows how this relation is modified once ex¬ 
tra gas physics is incorporated. We see that the increase in 
temperature associated with radiative cooling and/or pre¬ 
heating is precisely enough to bring the simulated relation 
into agreement with the observed one. 

In the Figures there are several clusters that lie well 
above the mean relation. These are mostly clusters for which 
there is significant velocity substructure; we indicate with 
open circles those clusters for which the mean gas and dark 
matter velocities within r25oo differ by more than 10 per cent 
of the rms velocity dispersion of the dark matter. 


5 DISCUSSION 

We have shown that simulations are capable of reproducing 
the observed relationship between mass and temperature in 
the inner regions of galaxy clusters. In particular, the mass- 
weighted temperature versus mass within a radius enclosing 
an overdensity of 2500 in our Radiative and Preheating sim¬ 
ulations agrees with the observed relation of ASFOl. 

There are a number of caveats, however. The temper¬ 
ature range of the simulations and the observations barely 
overlap; we have one cluster above 6keV, while ASFOl have 
only one below this temperature. Nevertheless, there is no 
reason to suppose that our results will not extend up to 
higher temperatures, though confirmation of this will have 
to await resimulations of clusters drawn from larger simula¬ 
tion boxes. 

Perhaps more pertinently, none of our simulations 
presents a fully realistic model of clusters, the Radiative 
model producing too much cooled gas and the Preheating 
model too little. However they both match the observed 
X-ray luminosity-temperature relation, because they both 
have a higher entropy within r25oo than does the Non- 
radiative simulation. This increase in entropy manifests it- 
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Figure 5. As for Figure^ but for the (a) Radiative, (b) Preheat¬ 
ing simulations. 


self as an inoroaoo in the tompomturo of the gao in the ir 


the clusters . One might expect , therefore , that re - 


parts o: 

alistic clusters that share the same entropy profile would 
predict the same temperature-mass relation. 

Unfortunately, the results presented in this paper and 
in ASFOl are of limited use to theorists who wish to predict 
the temperature function of clusters in order to constrain 
cosmology. This is because they need to relate the mass 
within the virial radius to the emission-weighted temper¬ 
ature of clusters. The prediction of masses at rsoo or larger 
radii from the X-ray observations is a harder problem than 
discussed here and will be investigated in a longer paper. 
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